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The  Sh- matrix  elements  derived  from  the  T-matrix  technique  allow  one  to  separate  the 
shape-dependent  parameters  from  size-  and  refractive-index-dependent  parameters. 
The  separation  also  allows  the  corresponding  surface  integrals  to  be  solved  analytically 
for  different  particle  shapes.  In  this  manuscript  we  give  analytical  expressions  for 
the  S/i-matrix  elements  that  contain  the  shape  information  for  cuboid-like  particles. 
We  find  very  good  agreement  between  the  results  obtained  using  the  Sh-matrix 
method  and  those  using  the  discrete-dipole  approximation  (DDA).  The  solution 
gives  significant  acceleration  of  calculations  dependently  on  the  size  parameter  of 
particles. 
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1.  Introduction 

Macroscopic  objects  and  small  particles  whose  shapes 
can  be  approximated  by  cuboids  (parallelepipeds  and 
prisms)  are  often  encountered  in  nature.  For  instance, 
many  different  crystals  can  be  described  as  cuboids  [1]. 
Thus,  the  problem  of  electromagnetic  wave  scattering  by 
such  particles  with  size  comparable  to  the  wavelength  is 
of  interest  in  different  fields  of  science  (e.g.,  astrophysics, 
atmosphere  and  ocean  optics,  biophysics)  and  engineering 
(e.g.,  radar  remote  sensing,  micro-technology). 

The  T-matrix  method  has  been  used  to  solve  the  light 
scattering  from  many  non-spherical  particles  [2-6]. 
It  is  especially  well  suited  for  atmospheric  aerosols  as 
orientation  averaging  can  be  performed  rapidly  using 
analytical  operations.  This  method  can  be  applied  to 
scattering  by  particles  of  arbitrary  shapes;  however, 
calculations  for  irregular  particles  are  difficult  and 
require  significant  computational  expense  to  compute 


*  Corresponding  author. 

E-mail  addresses:  petrov@astron.kharkov.ua  (D.  Petrov), 
gvideen@arl.army.mil,  gorden.videen@gmail.com  (G.  Videen). 

0022-4073/$ -see  front  matter  Published  by  Elsevier  Ltd. 
doi :  10.1 01 6/j.jqsrt.201 0.01. 024 


the  surface  integrals  necessary  to  find  the  elements  of  the 
T-matrix  [4].  Therefore,  development  of  approaches 
allowing  analytical  computations  of  these  integrals  is 
beneficial. 

Within  the  T-matrix  technique,  it  has  been  shown  that 
the  shape-dependent  factors  can  be  separated  from  the 
size-  and  refractive-index-dependent  factors  and  these 
are  contained  in  the  shape  matrix,  or  Sfr-matrix  [4,6-13]. 
The  elements  of  the  Sfr-matrix  are  expressed  in  terms  of 
surface  integrals.  Once  calculated,  the  T-matrix  of  parti¬ 
cles  that  have  other  sizes  and  refractive  indices  can  be 
found  through  analytical  operations.  In  many  cases  the 
surface  integrals  describing  the  Sfr-matrix  elements 
themselves  are  sufficiently  simplified  as  to  be  solvable 
analytically  (e.g.,  for  microlenses  [6],  Chebyshev  particles 
[7],  bi-spheres  and  capsules  [8],  merging  spheres  [9], 
finite  cylinders  [10],  corrugated  finite  cylinders  and 
capsules  [11],  cube-like  particles  [12],  and  ah  the  men¬ 
tioned  particles  with  a  concentric  layered  structure  [13]). 
In  this  manuscript,  we  present  an  analytical  solution  for 
the  Sfr-matrix  that  can  be  used  to  find  the  electromagnetic 
wave  scattering  from  cuboid-like  particles  with  smooth 
edges  and  vertices.  We  present  illustrations  of  the 
scattering  from  rectangular  parallelepipeds. 
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2.  S/i-matrices  for  right-cuboid-like  particles 

We  consider  a  family  of  particles  whose  shapes 
transition  from  ellipsoids  to  cuboids  expressed  as  follows: 

IWncM)- [(sin*)*  ((^)\ (^f)  +  (^ff * 

(1) 

where  6  and  </>  are  the  polar  and  azimuth  angles, 
respectively.  The  parameters  a,  b,  and  c  correspond  to 
the  half-lengths  of  the  particle  in  the  x,  y,  and  z  directions. 
When  a=b=c  the  particle  has  three-fold  symmetry  [12]. 
The  parameter  n0  characterizes  the  similarity  of  the 
particle  to  a  right  parallelepiped  (n0  should  be  even).  At 
n0= 2  the  particle  is  an  ellipsoid.  As  n0  increases,  the 
particle  becomes  less  like  an  ellipsoid  and  more  like  a 
parallelepiped.  Selected  particle  shapes  are  shown  in 
Fig.  1.  In  the  limiting  case  the  equation  for  a 
parallelepiped  is  as  follows: 

(2) 

The  Sh-matrix  elements  of  particles  whose  shape  is 
described  by  Eq.  (1)  can  be  expressed  analytically.  The 
derivation  is  extremely  long  and  tedious,  so  we  provide 
only  the  final  equations  describing  the  Sh-elements  for 
cuboid-like  particles  in  Appendix  A.  In  the  following 
sections  we  present  results  of  calculations  from  sample 
particles,  including  comparisons  with  DDA  simulations. 

3.  Calculations,  results  and  discussion 

In  our  calculations  all  particles  are  characterized  with 
the  refractive  index  m0  and  size  parameter  X=27ir/2, 


where  r  is  the  radius  of  the  circumscribing  sphere,  and  2  is 
the  wavelength.  Fig.  2  shows  maps  of  the  light-scattering 
intensity  distributions  in  the  forward  (a)  and  backward 
(b)  hemispheres  from  cuboid-like  particles  (m0=1.33  and 
size  parameters  Xa=12,  Xb=18,  and  Xc=24).  The 
distribution  of  the  degree  of  linear  polarization  is 
presented  in  Fig.  2c  (forward  hemisphere)  and  Fig.  2d 
(backward  hemisphere).  The  rows  in  these  figures 
correspond  to  particle  shapes  changing  from  ellipsoid 
(n0= 2)  at  the  top  to  more  cuboid-like  particles  (n0= 10)  at 
the  bottom.  The  maps  clearly  demonstrate  the  transition 
from  ellipsoid  to  cuboid.  As  the  particle  becomes  more 
like  a  cuboid,  the  ellipsoidally  shaped  scattering  fringes  in 
the  forward  direction  become  more  jagged  and  bright 
bands  appear  in  the  patterns  in  planes  perpendicular  to 
the  surface  facets.  Differences  are  enhanced  in  the 
backscattering  hemispheres.  The  columns  correspond  to 
the  fixed  orientation  of  particles  at  different  angles  cp  of 
electromagnetic  wave  incidence.  At  (p  =  0°  the  incident 
wave  propagates  parallel  to  the  normal  of  the  largest 
cuboid  face.  At  cp=45°  the  incident  wave  propagates  at  an 
angle  45°  to  this  normal,  and  at  (p= 90°,  the  incident  wave 
propagates  parallel  to  the  normal  to  the  smallest  cuboid 
face.  In  the  forward-scattering  direction  we  can  see  that 
the  fringe  spacing  is  inversely  proportional  to  the  size 
parameter  of  the  particle  cross-section,  with  larger  fringe 
spacing  occurring  in  the  plane  corresponding  to  the 
smaller  particle  dimension. 

Fig.  3  shows  the  normalized  intensity  and  degree  of 
linear  polarization  phase  functions  of  cuboid-like  particles 
averaged  over  orientation.  Curves  are  shown  for  particles 
corresponding  n0=2  (solid  lines),  n0=4  (dot  lines),  n0=6 
(dash-dot  lines)  and  no=10  (lines  with  points),  whose  size 
parameters  are  Xa=12,  Xb=18,  Xc=24  and  refractive  index 
is  m0=1.33.  As  can  be  seen,  the  parameter  n0  noticeably 


no  =  2  no  =  4  "«  =  6 


Fig.  1.  Examples  of  cuboid-like  particles  at  a=l,  b= 0.75,  c= 0.5,  having  different  n0. 
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Fig.  2.  Maps  of  the  light-scattering  intensity  distributions  of  cuboid-like  particles  (m0=1.33  and  size  parameters  Xa=  12,  Xb=  18,  Xc=24)  in  the  forward  (a) 
and  backward  (b)  hemispheres.  The  degree  of  linear  polarization  is  presented  in  (c)  (forward  hemisphere)  and  (d)  (backward  hemisphere).  The  rows  in 
these  figures  correspond  to  particle  shapes  changing  from  ellipsoid  (n0= 2)  at  the  top  to  more  cuboid-like  particles  (n0=10)  at  the  bottom.  The  columns 
correspond  to  the  fixed  orientation  of  particles  at  different  angles  cp  of  electromagnetic  wave  incidence.  At  cp=0°  the  incident  wave  propagates  parallel  to 
the  normal  of  the  largest  cuboid  face.  At  (p=45°  the  incident  wave  propagates  at  an  angle  45°  to  this  normal,  and  at  (p  =  90°,  the  incident  wave  propagates 
parallel  to  the  normal  to  the  smallest  cuboid  face. 


influences  the  phase  curves  of  intensity  and,  especially, 
the  polarization  degree.  In  particular,  at  small  phase 
angles  a  (a=7i-$)  the  ellipsoid  does  not  have  negative 
polarization  branch;  whereas,  one  is  present  in  all  cuboid¬ 
like  particles.  Fig.  4  presents  the  normalized  intensity  and 
degree  of  linear  polarization  of  cuboids  averaged  over 


orientation  when  one  of  the  three  dimensions  is  varied. 
Two  size  parameters  are  kept  constant  atXa=Xfa=32,  while 
the  third  size  parameter  is  varied  from  Xc=2  (plates)  to 
Xc=32  (cube-like).  The  relatively  large  size  parameters 
account  for  the  high-frequency  structure  in  the  phase 
curves.  While  all  particles  demonstrate  enhance  in  the 
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Fig.  3.  The  normalized  intensity  and  degree  of  linear  polarization  of 
cuboid-like  particles  averaged  over  orientation.  Curves  are  shown  for 
particles  having  n0= 2  (solid  lines),  n0=4  (dot  lines),  n0=6  (dash-dot 
lines)  and  no=10  (lines  with  points),  whose  size  parameters  are  Xa=12, 
Xfa=18,  Xc=24  and  the  refractive  index  is  m0=1.33. 

backscattering  intensity  and  a  significant  negative 
polarization  branch,  a  significant  positive  polarization 
branch  is  only  present  when  the  third  size  parameter  is 
small.  Fig.  5  shows  the  normalized  intensity  and  degree  of 
linear  polarization  for  crystals  of  two  minerals,  olivine  and 
topaz,  which  are  averaged  over  orientations.  For  olivine, 
the  ratio  of  the  dimensions  are  a:b:c=0.466:l:0.589  and 
m0=1.67;  for  mineral  topaz  a:b:c=0.528:l :0.954  and 
m0=1.63.  The  size  parameters  of  these  crystals  were 
chosen  to  be  equal  to  those  of  the  sphere  of  equal  volume, 
having  X=5.0.  While  the  intensity  phase  functions  do  not 
vary  significantly  between  the  two  types  of  particles,  the 
polarization  phase  functions  are  significantly  different. 
These  calculations  can  be  used  as  a  benchmark  in  optical 
studies  of  the  minerals. 

To  verify  our  theory  and  code,  we  have  compared  Sh- 
matrix  calculations  with  those  made  using  the  discrete- 
dipole  approximation  (DDA)  [14].  Fig.  6  shows  the 
dependence  of  intensity  and  degree  of  linear 
polarization  as  a  function  of  the  scattering  angle  $  for 
cuboid-like  particle  with  the  following  parameters: 
n0= 20,  Xa=5.0,  Xb=7.5,  Xc=10.0  and  m0=1.5  +  0.1i.  The 
particles  are  oriented  so  that  the  incident  electromagnetic 
is  perpendicular  to  the  largest  cuboid  face.  Points  and 
solid  lines  correspond  to  the  calculations  carried  out  by 


Fig.  4.  The  normalized  intensity  and  degree  of  linear  polarization  of 
cuboids  (square  plates)  averaged  over  orientation,  whose  basis  is  a  square 
with  Xa=Xb=32  and  Xc  that  varies  from  2  to  32;  m0=  1.33  and  no=10. 


Fig.  5.  The  normalized  intensity  and  degree  of  linear  polarization  of  crystals 
of  olivine  and  topaz,  averaged  over  orientations.  For  olivine  a:b:c=0.466: 
1:0.589  and  m0=1.67.  For  topaz  a:b:c=0.528:l:0.954  and  m0=1.63.  The 
volume  of  these  crystals  was  chosen  to  be  equal  to  X=5.0  spheres. 
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Fig.  6.  The  intensity  and  degree  of  linear  polarization  of  a  cuboid-like 
particle  in  a  fixed  orientation  illuminated  by  an  incident  electromagnetic 
wave  propagating  perpendicular  to  the  largest  cuboid  face.  The  following 
parameters  are  used:  n0= 20,  Xa=5.0,  Xb=7.5,  Xc=10.0  and  mo=1.5  +  0.1i. 
Points  correspond  to  the  calculations  carried  out  with  the  DDA  method 
in  Zubko’s  version  [14]  and  solid  lines  are  calculations  with  the  Sh- 
matrices  method. 


means  of  DDA  and  5/i-matrices  methods,  respectively.  As 
can  be  seen,  the  curves  coincide  very  well. 

It  also  should  be  noted  that  the  analytical  solutions 
give  significant  acceleration  of  calculations  as  compared 
to  the  direct  computation  of  the  surface  integrals 
dependency  on  the  size  parameter  of  particles;  however, 
the  primary  time-consuming  portion  of  the  calculations 
remains  computing  the  inverse  matrix  that  is  common  for 
all  variations  of  the  T-matrix  method. 


4.  Conclusion 

We  have  used  the  5h-matrix  approach  to  derive 
analytical  solutions  for  a  family  of  shapes  representing 
cuboid-like  particles.  Ellipsoids  and  pure  cuboids  are  two 
extreme  cases  from  this  family.  For  cuboid-like  particles 
edge  diffraction  bands  are  a  dominant  characteristic  of  the 
scattering  patterns  giving  cross-like  structure  in  the 
angular  distributions  of  the  intensity  and  degree  of 
polarization.  Although  analytical  solutions  significantly 
accelerate  calculations  of  scattering  characteristic,  the 
main  time-consuming  step  remains  computing  the  in¬ 
verse  matrix. 
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Appendix  A 

In  this  appendix  we  present  general  expressions  for  the  Sh-matrix  elements  of  the  T-matrix 
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(—Ijm'-m  +  n-l  +  /< 

Shmnmn',k  = 22k  +  „--„  +  i  Ann' 


\m'-m-n+ 1  +  k 


-n^Ln'(2  k-n+n') 


,“mnm'n'vt"v  11  1  / 

+  ^ mrnn'n '  (2/C-H  +  717  - 1 )  +  L^nm'n'  (2/C-H  +  n7) 


Shmnm'n',k  =  ^2k+n'-n  +  2  Kn'l^nm' +  2), 


5^211  -  fTT  (~1)?T .  ^ 

3 nmnm'n',k ~  171  ^ 2k+n-n  +  t 

\m' -m-n-\  +k 


(n'  +  l)/®m,„,(2k-n+n') 

+C°m'n'(2fe— n  +  n'-l  )-iL("ln,(2k-n  +  n') 


/ _ -1  \m  —m—n—  i  +  k 

,v  Am'I%lm,„,(2k-n+n' +  2), 


C/j212  -i7r{~') _ 

,lmnm'n',k  ^k  +  n'-n  +  l 


0,221  _  Azll _ A 

3 nmnm'n',k-n  ^ 2/<+n'-n+l 


r-n(n'  +  l)/^m,„,(2k-n+n'-l) 


+  (n7  +  l)^mnm' 


-nL(U{ 

u^m'n'mn 


c  1,222  _  4 

^ nmnm'n',k  ~  ~n  ^k+rt-n+2  ^ 


cu223  _  —  v  A1  4 

*nmnm'n',k  ~  71  2k  +  n'-n  +  2  ^ 


-mnm'n'(2k— n- +  n'-2) + i(nr  + 1  )L(n5J(n,n,  (2k-n  +  n'  - 
<00)  (2k-n  +  n'-2) + inLi'JJ,mn(2k-n  +  n'-l) 

-nCm'n(2fe-n  +  n'  +  l) 

+ Cw  (2fe-  n  +  n')  +  iL'„],,ln,n.  (2k-n+ri  +  l) 

(>T  +  +Cnm,n.(2k-n+n+l) 

+Cmn(2(f-n+ n')-iCl>m„(2k-n  +  n'  + 1) 


/  -i  \m'-m-n  + 1  +  k 

Shfnnnvn'.k  =  71 - Z^Fwr^n  +  3 - Ann'Im)„mni(2k-n  +  ri  +3). 


The 


nmnm'n',k~  ^ 2/<  +  n/-n  +  3  nm'Jmnm'n^z,,v 

following  designations  are  used  in  formulae  (A.1)-(A.18) 


CmvCZ)  =  ™ 


n'y/(n'  +  l)2-|T1'2  ^9)  (n'  +  l)Vn'2-m'2  .(g) 

2n'  +  l  1mnm'n'  +  1^  2n'  +  l  1mnm'n'~  lv6/ 


+  m7 


n\An  +  l)2-m2  9)  (n  +  l)Vn2-m2  (9) 

2n_j_4  1mn+\m'n''A)  2n  +  l  1mn-\m'n,'A) 


(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A.  10) 

(All) 

(A.  12) 

(A.  13) 

(A.  14) 

(A.  15) 

(A.  16) 

(A.  17) 

(A.  18) 


(A.  19) 


Cm'n'(z)  =  |mm'lCm,„,(z)+  (2„,  +  1)(2n  +  i)  [(n'  +  1)v/ n'2-m'2(n  +  1)V n2-m2C-iwn'-i(z) 

+  n'  (n'  +  l)2-m'2n \J{n  +  \f-m2f®n + ,  m,n,  + ,  (z)-(n'  + 1 )  Vn'2-m'2n  f{n  +  \f-m2t^ + ,  (z) 
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-ri  \J  (n'  +  l)2-m'2(n  + 1  + ,  (z) 

where 


£(0°)  ("2)  —  mrlr(0)  (z) 

Lmnm'n'W  —  mi  xnnva'n'^'-’ 


r  (01) 


,,n'(Z)_  2n'  +  l 


nV(n'  +  l)2-m'2,'» + 1  (z)-(n'  +  (z) 


Lmnl'n'(Z)  =  m'I(£m,n,(Z), 

Lrnl- * (z)  =  2^  \/(n'  + 1  )2_m'2Cm'n' + , (z)-(n'  + 1 ) Vn'2-m'2CmW-, (P 


and 


(A.  20) 

(A.21) 
(A.22) 
(A.23) 
(A.  2  4) 


n-|m|  _ 

. ,-ir"'2„2„„!V(»-»l)!(n  +  |ml,!»V(»'-ln.'l)!(n'  +  ln.'l)!  E  tK„-t|.(„-M-t)K|mH.|t)! 


(-1)* 


n'-|m'|  /i-vK1 

x  V  7— - =— — - — - =— [/(C1)(m'-in,2n-2/f-|m|+2n'-2k'-|m'|-l,2k 

^  k'!(n'-k')!(n'-|m'|-k')!(|m'|+k')!  1  1 

+  |m|+2k'  +  |m'|-l,z)], 

n-|m|  /_  ixk 

'  ™"-(Z)  =  (— l)n+n  Sm£m<n!^(n-|m|)!(n  +  |m|)!n'!-y/(n'— |m'|)!(n'  +  |m'|)!  g  k!(n_k)!(n_|m|_k)!(|m|+k)! 


(A.25) 


n'-|m'| 


k  =  0 


V  -p— - — — - — - — - p— [/(C2)(m/-m,2n-2/c-|m|  +  2n/-2/c/-|m/|,2/<:+|m|  +  2/c/  +  |m/|,z)] 

(A.  2  6) 


fc'  =o 


_ n-|m|  _ 

Cm,n,(z)  =  (-l)n+n 'SmSm-n!^/(n-|m|)!(n+|m|)!n'!^/(n'-|m'|)!(n'  +  |m'|)!  g^  k!(n_k)!(n_t~y_k)!(|m|+fc). 


(-If 


E 

K  =0 


(-If 


where 


fa\m,r],v,z)=  £  "°fc  ^E(-lf'C^(m,p') 

p  =  0  Pm  p'  =  0 

p-p'  Ql<  n0l< 

x  E  T^f  E  (-2fC%okQ(v  +  n0p'  +  2k',ti+n0p'X 

k  =  0  L  k'  =  0 

v.z) -2  f  (^  +  1)(^)  r(^-<"2)  £ 

p  =  0  P-  p’  =  0 


(A.  2  7) 


(A.  2  8) 


p-p'  n0k  +  n0-1 

E  (-2f'CQ(v  +  n0p'  +  2k'  +  l,;7+nop'  +  l)-2"o-1 

/<  =  0  k'  =  0 

~  ((z+l/n0)+l)(z+l/n0)...  ((z+l/n0)-p+2)  p  ,  , 

X  E- - - - s - X - J-'£(-\ycr<P(m,p'  +  1) 

p  =  0  P'  p'  =  0 

p-p'  £k  n0k+ 1 

E^nSr  E  (-2f'C^fl(v+n0p'+2k'+no-l,p  +  n0p'  +  no-l), 


cn0/<  + 

k  =  0  k  =  0 

oo  fz+l_n  +  2^|  P 

fiC3\m,rj,v,z)  =  -2n°^^ - )\n0j---yn0  - Z  £(_1  fCP0(m,p') 

p  =  0  p'  =  0 


p-p'  (j/c  n0/c  +  n0-l 

E  rn04~n0-i  E  (-2f'C^(v  +  nop'  +  2k'  +  no,p  +  n0p'  +  no) 


where 


rn0k  +  n0 

k=  0  L  k  =  0 


(A.  2  9) 


(A.  30) 


p  mexp[nt(m-n0p+nok-l/2)](  l+exp[ntm](-l)"oP) 
<I>(m,p)  =  £  c£ - ^ - 7 


fc  =  o 


Qn„kf,noP-nolf  2noP(n0k+  1  )B  +  1 ,  +2) 


2Fl(nofe-„oP,!!^P;l+2^°P  +  m;-1 


(A.31) 
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Q(f],  V)  = 


2r(^+i) 


0(m,p)=  Y,C\ 


inexp[in(m-n0p + nok-2/2)]  (l  +  exp[inm](- 1  )"oP) 


kZ o  P  ano/<+»o-i^nop-no/<+i2noP+no(n0/<:+no)B^(noP  +  no-m/2)  +  l,(2no/c-noP  +  no  +  m/2)^ 

x2Ft  (ji0k-n0p-l,  (m-riop-v/2);  (2n0k-n0p  +  n0  +  m/2);  -l) 

p  ^  inexp[in(rn-n0p+n0k-n0/2)](^\  +  exp[mm](-l)noP) 

k5)  P  anofc+1bnoP-nofc+no-i2noP+no(no/c+2)B^(n0p  +  no-m/2)  +  l,(2no/c-nop-no-m  /2)+2) 
x2Fi  (tiok-nop-no  + 1 ,  (m-n0p-n0/2);  (2n0k-n0p-n0-m/2)+2; -l) 
where  2Fi  is  the  Gaussian  hypergeometrical  function 


2F1(a1,a2;b1;z)  = 


Hb,) 


ria,  )r(a2)Z-(l 

B(n,m )  is  the  Beta-function 


F  (a,  -fnl/Xaz-nOz11 3 4 5 6 
r(bi+n)  n!’ 


B(n,  m)  = 


mnm) 

r(n+m )  ’ 


(A.32) 


(A.33) 


(A.34) 


(A.35) 


J  1,  m>  0 

m<  0, 


(A.36) 


Ann'  = 


(2n'  + 1)  /  (2n  +  l) 


47HT (n7  + 1 )  y  47in(n  + 1 ) 


and  Cnm=(n!/m!(n-m)!)  are  the  binomial  coefficients. 


(A.37) 
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